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CORRECTIONS  FOR  FREQUENCY  RESPONSE  DIFFERENCES 
IN  MULTIPLE  RECEIVER  ARRAYS 


by 

J.L.  Robinson  and  R.W.  Jenkins 


ABSTRACT 

' ''  In  interference-cancelling  antenna  array  systems  where  the 
signals  are  combined  after  passing  through  receivers,  the  differential 
frequency  responses  of  the  receivers  act  to  limit  the  null  depth  and 
therefore  the  degree  of  interference  rejection  attainable  for  finite- 
bandwidth  interference.  An  experimental  four-element  HF  antenna 
array  is  used  as  an  illustration;  the  measured  differential  frequency- 
dependent  phase  and  amplitude  responses  of  the  receivers  are  pre¬ 
sented  and  the  corresponding  limits  to  null-depth  derived,  for  a  3 
KHz  bandwidth.  Correction  techniques,  including  an  exact  software 
procedure  and  a  hardware-implementable  approach  to  linear  phase 
correction  are  presented,  and  their  requirements  and  anticipated 
performance  considered 


1 .  INTRODUCTION 

Interference-cancelling  antenna  arrays  operate  by  placing  nulls  in 
the  direction  of  unwanted  signals  while  at  the  same  time  maintaining  non-zero 
gain  in  other  directions.  This  is  done  by  weighting  i.e.  adjusting  the 
amplitude  and  phase  of  the  signals  received  from  the  various  array  elements 
and  then  combining  them.  Unwanted  signals  are  adjusted  so  as  to  cancel  upon 
combination. 

In  many  arrays,  the  signals  pass  through  receivers  before  they  are 
weighted  and  combined.  In  such  systems,  differences  in  the  frequency  responses 
of  the  receivers  cause  a  frequency  variation  of  the  relative  signal  amplitudes 
and  phases  between  receivers.  Frequency-independent  weighting  which  allows  a 


signal  component  at  one  frequency  to  completely  cancel  itself  upon  combination 
will  therefore  not  be  effective  for  signal  components  at  other  frequencies, 
and  a  limit  to  the  null  depth  obtainable  for  non  zero-bandwidth  signals  is  set. 


An  estimate  for  this  null  depth  is  derived  in  the  Appendix.  The 

result  is 


p  /p  = 
N  A 


+  a. 


(1) 


where  PN  is  the  array  gain  in  the  direction  of  the  interfering  signal,  PA  is 
the  average  array  gain,  and  Op2  and  o^2  are  the  variances  in  receiver  gain 
(expressed  as  a  dimensionless  ratio)  and  phase  delay  (in  radians2),  taken 
relative  to  the  mean  frequency  response  of  the  receivers.  Some  representative 
null  depths  for  various  values  of  Op  and  0^  are  given  in  Table  1.  These 
values  can  also  be  considered  as  the  limits  to  Op  and  o^  that  must  not  be 
exceeded  in  order  to  achieve  a  certain  level  of  interference  cancellation. 


TABLE  1 


Achievable  Null  Depths,  for  Various  Equal-Contribution  Amplitude  and  Phase  Variances 


PN/PA  <dB> 

°P 

0  p  (degrees) 

-20 

.071 

~ i 

4.1  I 

-30 

.022 

1.3 

-40 

.007 

0.4 

-50 

.002 

0.1 

In  this  paper  we  consider  the  magnitude  of  the  above  described 
differential  frequency  response  effect  and  techniques  for  reducing  it.  The 
case  of  an  experimental  four-element  HF  array  using  good  quality  commercially 
available  receivers  and  employing  adaptive  interference  cancelling  at  base¬ 
band  is  used  as  an  illustration.  A  series  of  measurements  of  the  amplitude 
and  phase  response  of  the  receivers  in  the  experimental  array  is  presented 
and  the  null  depths  achievable  for  broadband  Interference  are  calculated.  Two 
correction  techniques  are  described:  a  precise  software-implementable 
correction  technique  that  uses  Fourier  transforms,  and  an  approximate 
correction  technique  involving  a  linear  phase  correction  using  time  delays. 
Required  processing  times,  expected  corrections  and  some  general  conclusions 
are  discussed. 


2.  MEASUREMENTS 


2.1  DESCRIPTION 

The  equipment  configuration  for  the  measurements  is  shown  in  Figure  1 
The  input  RF  signal  was  generated  by  a  frequency  synthesizer  (Hewlett  Packard 
3335A),  suitably  attenuated,  split  four  ways,  and  routed  through  the  four 
receivers.  Each  receiver  (RACAL  6790)  was  tuned  to  the  input  RF  and  set  for  a 
6  kHz  bandwidth.  The  receivers  were  driven  by  a  common  reference  signal,  so 
that  the  relative  phases  of  the  input  signals  were  maintained.  The  455  kHz  IF 
outputs  were  mixed  down  to  baseband  and  split  into  in-phase  (I)  and  quadrature 
(Q)  channels  in  the  hybrid  quadrature  mixers.  The  455  kHz  mixing  signal  was 
provided  externally  by  a  frequency  synthesizer  (HP  3335A) .  The  eight  output 
signals  were  routed  through  low  pass  filters  with  an  upper  cut-off  frequency 
of  1.5  kHz  (Precision  Instruments  filter  unit).  One  line  was  chosen  as  a 
reference  and  connected  via  the  channel  A  input  to  the  spectrum  analyzer 
(HP  3582A).  The  remaining  lines  were  connected  in  turn  to  the  channel  B  input 
The  spectrum  analyzer  provided  a  digital  readout  of  the  amplitude  for  all 
eight  lines  and  relative  phase  measurements  of  the  last  seven  lines  with 
respect  to  the  first  line.  The  receivers  were  tuned  to  a  center  frequency  of 
11  MHz  and  the  input  signal  was  stepped  in  100  Hz  increments  from  10.9985  MHz 
to  11.0015  MHz. 


LOW 

QUADRATURE  PASS 
RECEIVER  /MIXER  FILTER 


Figure  1.  Measurement  Equipment  Configuration 


2.2  RESULTS 


Figure  2  shows  the  amplitude  response  of  the  eight  output  channels 
across  the  frequency  band.  There  is  a  difference  in  the  absolute  magnitude 
of  the  receiver  outputs  of  at  most  4.5  dB,  and  a  much  smaller  variation  with 
frequency  in  the  ratios  of  amplitudes  over  the  receiver  bandwidth  (at  most 
0.2  dB) .  This  latter  feature  is  seen  clearly  in  the  bottom  of  Figure  2.  In 
this  figure  the  distribution  about  the  mean  of  the  amplitude  of  channel  3 
relative  to  the  amplitude  of  channel  1  is  shown.  The  variance  for  these 
results  has  been  calculated  to  be  cfp2  =  3.36x10”^.  The  total  variance  for  all 
channels  relative  to  the  average  channel  amplitude  response  was  also  calcu¬ 
lated,  giving  Op2  =  1.051x10”^. 


FREQUENCY  (KMi) 

Figure  2.  Amplitude  Responses  of  the  Eight  Outputs  and  the  Ratio  of  Amplitudes  for  Output  3  Relative  to  Output  1. 


Figures  3(a),  (b) ,  (c)  and  (d)  show  the  relative  phase  response  of 
channels  2  through  8  with  respect  to  channel  1.  Figure  3(a)  shows  the 
measured  phase  difference  between  the  I  and  Q  outputs  of  receiver  1,  which  as 
expected,  is  approximately  90  degrees  and  shows  little  variation  across  the 
pass  band.  Figures  3(b)  -  (d)  show  a  characteristic  monotonic  decrease  in 
relative  phase  across  the  ±1.5KHz  frequency  band.  Least  mean  square  straight 
line  fits  to  the  data  have  been  calculated  and  are  included  in  these  figures. 
The  fluctuations  of  the  data  about  these  straight  line  fits  are  approximately 
the  same  for  the  I  and  Q  outputs  from  the  same  receiver.  The  straight  line 
fits  are  further  discussed  in  Section  3.2.  The  total  variance  of  all  phase 
data  points  with  respect  to  the  average  frequency  response  was  estimated  from 
the  phase  data  to  be  o^2  =  5.3x10“^  (radians^). 


FREQUENCY  UMi) 


Figure  3.  Relative  phase  response  of  the  outputs  with  respect  to  output  1.  M  is  the  slope  and  B  is  the  zero 
intercept  of  the  least  mean  square  straight  line  fit  to  the  data. 
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It  follows  from  Equation  (1) ,  using  the  value  of  Op2  and  CU2  given 
,  that  the  differential  frequency  responses  of  the  receivers  will  limit 


above,  that  th 
null  depths  to 


PN/PA  =  °p2  +  °<p2 

=  l.OSlxlO-4  +  5.3xl0~3  =  5.4xl0“3,  or  -22.7dB 
with  the  largest  contributing  factor  to  the  limit  being  the  phase  differences. 

3.  CORRECTION  TECHNIQUES 

Two  techniques  are  described  for  making  corrections  to  the  differ¬ 
ential  frequency  response  to  assure  satisfactory  null  depths.  First,  a 
completely  general  computational  approach  that  uses  frequency  domain  amplitude 
and  phase  measurements  is  described  in  detail.  This  approach  can  be  used  if 
sufficient  signal  processing  power  is  available.  Second,  a  hardware  approach 
is  described  that  staggers  the  sampling  times  and  can  be  used  when  a  linear 
correction  to  the  phase  error  is  sufficiently  accurate. 

3.1  FOURIER  TRANSFORM  TECHNIQUE 

A  correction  for  the  difference  in  the  frequency  response  of  one  or 
more  receivers  relative  to  a  standard  receiver  can  be  made  by  directly 
correcting  the  data  in  the  frequency  domain.  This  technique  requires  a 
processor  that  has  a  real  time  Fourier  transform  capability.  A  finite  sequence 
digital  compensation  filter  is  constructed  using  frequency  response  measure¬ 
ments  of  the  receivers  and  is  kept  available  in  the  signal  processor  memory. 

The  correction  is  accomplished  by  a  discrete  linear  convolution  of  the  finite 
filter  sequence  with  the  indefinite  length  incoming  signal  sequence  [1].  This 
is  best  implemented  by  a  multiplication  in  the  frequency  domain.  The  entire 
correction  procedure  is  accomplished  in  software  and  should  be  designed  to  be 
as  fully  automated  as  possible  for  operator  convenience. 

The  construction  of  the  frequency  domain  sequence  representing  the 
compensation  filter  is  based  on  calibration  data  obtained  using  the  experi¬ 
mental  configuration  given  in  Figure  1  and  following  the  measurement 
procedure  given  below: 

1.  Set  the  synthesizer  to  an  appropriate  starting  frequency 

(fo>- 

2.  Distribute  the  sine  wave  signal  from  the  synthesizer 
equally  to  all  receiver  inputs. 

3.  Allow  sufficient  time  for  steady  state  conditions  to 
become  established . 


v  .V  ,s  \  ,v  - 

^  v  v  *. 


-1  . 


.*■  .V V *V-\' V-\*  Vw\v. 

v  v  ■.  v  v  »•  v  'a'.*  v  v  v  * 
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4.  Perform  A/D  sampling  and  Fourier  transformation  of 
receiver  output  to  obtain  amplitude  and  ; base  measure¬ 
ments  at  fQ. 

5.  Repeat  steps  (l)-(4)  for  all  frequencies  across  the 
response  band  of  the  receivers,  in  steps  of  width  Af. 

In  accordance  with  step  if 4,  the  digitized  receiver  outputs  give  a 
time  sequence  x^t^) ;  i=l,...,N  at  each  input  frequency  for  each  output 
channel  n.  A  discrete  Fourier  Transform  (DFT)  is  computed  at  each  input 
frequency  with  the  N  samples  from  each  receiver  output,  thereby  providing  the 
complex  frequency  domain  representation  of  the  output.  At  an  arbitrary 
frequency  fk  for  example. 


xn(fk>  =  1[1*n(ti)a“'(--12rtkti)  ‘  An(fk)exp('J<(,n(fk)) 


where  x^ffc)  is  the  output  from  channel  n  for  the  input  signal  at  frequency 
fk  and  ^(f^)  and  4>n(fk)  are  the  amplitude  and  phase  respectively  of  the 
frequency  domain  representation  of  the  channel  n  output  signal. 

With  one  output,  m,  chosen  as  the  standard,  the  filter  coefficient 
at  frequency  fk  for  any  other  output  n,  is  defined  by 

Vfk>  (-KVV-W) 

where  and  4>n,-<t>n  are  the  relative  amplitude  ratio  and  phase  difference 

of  the  nth  output  with  respect  to  the  mth  output.  It  is  essential  that  the 
sampling  rate  used  during  calibration  is  matched  to  the  sampling  rate  of  the 
Implementation.  This  is  ensured  by  the  following  procedure  (the  requirement 
for  this  matching  procedure  is  discussed  below) : 

1.  Begin  with  Kn(fk)  obtained  during  the  calibration  at 
Q  points  with  frequency  spacing  Af. 

2.  Expand  to  M**2Q-1  points  using  the  symmetry  relation: 

Kn(fk)  -  MMAf-flc)  • 

3.  Compute  an  M  point  inverse  FFT;  the  Impulse  response 
has  length  M  and  spacing  At  •  1/MAf. 

4.  Insert  M  zeroes  onto  the  end  of  the  impulse  response. 

5.  Compute  a  2M  point  FFT;  the  new  K^f)  has  Mfl  unique 
points  with  frequency  spacing  Af'  -  l/2MAt  ■  Af/2. 


.*.  •  +  '  **  »  •  Ajis  *  *  V  •  A  >N  .V  ►  •  .  •  *  -  >  .  ■  .  V  •  *  «  .  *  ^  A  .  •  ,*  .*»  .  *  A  .S'  ^  . 
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This  filter  coefficient  is  now  the  desired  factor  and  the  multiplication 


WW  ■  w 


makes  the  desired  correction,  where  K^f^)  is  stored  in  the  processor  memory 
and  XjjCffc)  is  the  DFT  component  at  frequency  f^  of  the  signal  from  receiver 
output  n.  This  will  correct  the  frequency  response  error  relative  to  channel 
m. 


The  discrete  linear  convolution  procedure  for  correcting  variations 
in  the  frequency  response  of  the  receivers  involves  the  computation  of  DFTs 
on  consecutive  blocks  of  the  incoming  sampled  time  signal.  These  blocks  are 
multiplied  by  K^f^)  in  the  frequency  domain  and  then  an  inverse  DFT  is 
computed  giving  a  finite  output  sequence.  The  output  blocks  must  be 
carefully  assembled  to  give  the  correct  indefinite  length  output  time 
sequence.  One  iteration  of  this  operation  as  it  would  be  applied  to  each 
channel  using  the  overlap  and  add  procedure  [2]  involves  the  following  steps: 

1.  The  input  of  M  real  time  samples  is  augmented  with 
M  zeroes. 

2.  A  2M  point  FFT  of  the  time  sequence  is  computed:  this 
gives  for  output  n,  X^Cf)  with  M+l  unique  points. 

3.  Compute  X^Cf^)  =  ^(f^K^f^)  (for  all  M+l  frequency 
points;  generate  another  M-l  points  using  symmetry). 

4.  A  2M  point  inverse  FFT  of  X'n(f)  is  computed  to  get  2M 
time  domain  points. 

5.  The  first  M  points  of  the  time  domain  sequence  are  added 
to  the  last  M  points  from  the  previous  iteration;  the 
second  M  points  are  saved. 

6.  Repeat  (1)  -  (5)  with  the  next  M  time  samples. 

Care  most  be  taken  in  choosing  the  calibration  frequency  step  size 
and  the  time  sampling  rate  to  avoid  seriously  distorted  results  due  to 
spectral  leakage.  The  step  size  in  the  filtering  operation  must  be  compatible 
with  the  step  size  in  the  calibration  measurement  [2].  This  is  ensured  by 
the  choice  of  compatible  sampling  rates  and  a  2ero  padding  procedure  that  is 
used  in  contructing  the  filter  K^f^).  When  the  time  sampling  rate  is  1/At 
and  the  Incoming  time  signal  is  being  processed  in  groups  of  M  samples  the 
frequency  spacing  in  the  calibration  measurement  is  Af  =  1/MAt. 

Most  of  the  processing  time  required  by  the  correction  algorithm  is 
taken  up  with  Fourier  transforms  and  multiplication.  The  implementation  of 
the  correction  for  one  channel  at  a  sampling  rate  of  5K  samples/second  and 
block  size  of  M*1024  samples,  would  require  approximately  5  iterations  per 
second.  For  eight  channels  this  Increases  to  40  iterations  per  second  or  a 
time  of  25  milliseconds  per  convolution.  Tests  [3]  of  the  processing  speed 
of  one  array  processor  (AP  400)  have  found  that  a  1024  point  complex 


convolution  operation  (FFT+multiplication+IFFT)  can  be  computed  within  a  time 
of  20  milliseconds.  Therefore  it  is  concluded  that  a  differential  frequency 
response  correction  based  on  the  Fourier  convolution  filtering  technique  can 
be  implemented  using  available  state-of-the-art  high  speed  processors. 


3.2  LINEAR  CORRECTION 

The  results  presented  in  this  paper  for  a  four-element  array  indicate 
that  the  extent  to  which  signals  can  be  nulled  is  limited  more  by  variations 
in  the  phase  response  of  the  channels  than  variations  in  the  amplitude 
response.  Also  the  overall  variation  of  phase  as  seen  in  Figure  3,  is  largely 
a  linear  function  of  frequency,  at  least  over  the  3  kHz  bandwidths  considered. 

A  differential  phase  response  that  varies  linearly  with  frequency 
can  be  considered  as  a  relative  time  delay  between  the  receiver  output  lines. 
Suppose  a  signal  exp(j2rft)  is  incident  on  the  inputs  to  lines  1  and  n,  and 
the  effect  of  the  two  lines  is  to  introduce  a  time  delay  T^,  (i=l,n),  plus  an 
arbitrary  amplitude  gain  and  phase  change.  The  output  signals  are  then: 


A^exp(  j  (2Trf  (t— T^)+0^)  )  for  line  1 
A^exp( j  (2rrf  (t-T^J+S^)  )  for  line  n 
The  differential  phase-response  is  then  given  by: 


Am*  -  -2irf(T1-Tn)+ei-0n  =  2ir(AlnT)f+(e1-en) 

where  A^nT  =  Tj.  -  Tn  is  the  time  between  lines.  Relating  this  to  the  observed 
linear  dependence  A^n*  =  M^nf  +  B^n  where  M^n  is  the  slope  and  B^n  is  the 
zero  intercept,  we  find  the  differential  time  delay  to  be  A^nT  =  -Min/2ir. 

Table  2  gives  the  differences  in  time  delay  between  channels  2  to  8 
and  reference  channel  1,  derived  from  the  straight-line  fits  to  the  data  in 
Figure  3. 


The  maximum  observed  time  delay  is  of  the  order  of  20  microseconds, 
which  is  considerably  less  than  the  333  microsecond  Nyquist  complex  sampling 
interval  required  for  3  kHz  bandwidth  signals.  For  a  digital  receiving 
system,  time  delays  of  that  magnitude  can  be  readily  corrected  by  making 
small  adjustments  to  the  sampling  time  for  the  analog-to-digital  converter  at 
each  receiver  output.  A  possible  hardware  implementation  is  shown  in  Figure  4. 

Once  a  time  delay  correction  is  made,  the  linearly  dependent  portion 
of  the  differential  frequency  phase-response  will  be  removed.  From  the 
residual  differential  frequency  responses,  a  new  upper  bound  to  the  inter¬ 
ference  cancellation  achievable  can  be  found,  by  using  Equation  1.  We  get, 
for  the  measured  receivers  after  time  delay  correction 


a  2  =  2.52x10 
P 

?  -4 

and  o.  =  1.05x10  as  before,  so  that 

w 

-4 

P„/P.  =  3.57x10  ,  or  -34.5dB.  This  is  a  considerable  improvement  over 

N  A 

the  previous  limit  of  -22.7dB  found  for  the  uncorrected  measurements  in 
Section  2.2. 


TABLE  2 

Values  M  i  n  of  the  Slope  of  the  Straight  Line  Fits  to  the  Differential  Phase  Response  Results  of  Figure  3, 
and  the  Corresponding  Implied  Differential  Time  Delays  nT 


i 

| 


I 


Channel  # 

M1n 

AmT  ! 

1 

1 

0°/KHz 

0  microseconds  j 

2 

-0.25 

0.69  ! 

3 

-4.76 

13.22  j 

4 

-4.59 

12.75  I 

5 

-7.02 

19.50  | 

6 

-5.32 

14.78  i 

7  ! 

-2.79 

7.75  ; 

8 

-5.1 1 

14.19 

1 

etc 


Figure  4.  Possible  Hardware  Implementation  for  Linear  Correction  by  Staggered  Sampling  Times 


4.  DISCUSSION  AND  CONCLUSIONS 


Of 


Achievable  null  depths  for  an  operational  adaptive  antenna  are 
governed,  according  to  Equation  1,  by  the  differential  frequency  response  of 
Che  receiving  equipment  located  ahead  of  the  signal- summing  apparatus.  This 
equipment  includes  the  antennas  and  any  feed  lines,  preamplifiers,  filters 
and  receivers  which  precede  the  circuit  that  does  the  summation. 

The  measurements  reported  herein,  on  a  set  of  four  identical  good 
commercial  grade  receivers,  showed  that  without  correction  the  null  depth 
for  a  3  kHz  bandwidth  signal  in  an  array  using  these  receivers  would  be 
limited  to  22.7  dB.  It  was  observed  that  this  limit  was  mostly  due  to  the 
differences  in  phase  response  between  receivers  (17dB  above  the  amplitude 
response  limitation).  Further,  the  phase  differences  in  frequency  response 
were  mainly  linear  functions  of  frequency,  which  permitted  their  interpre¬ 
tation  as  frequency-independent  differences  in  receiver  time-delay.  If  these 
time  delays  were  corrected  for,  the  residual  phase-response  contributions 
would  be  about  the  same  order  of  magnitude  as  the  differential  amplitude 
response  contribution.  The  resulting  null  depth  would  be  34.5  dB  for  a  3  kHz 
bandwidth  signal.  A  linear  phase  difference  correction  is  easily  Implemented 
in  a  digitally  sampled  system  by  adjusting  the  relative  sampling  times. 

Another  factor  which  limits  the  attainable  null  depth  is  the  siting 
of  antenna  elements.  A  previous  set  of  measurements  [4]  in  one  instance 
indicated  worst-case  null  depths  due  to  antenna  element  siting  limited  to' 

35  dB.  Therefore,  the  presently  considered  receivers,  together  with  a  time 
delay  correction,  would  not  greatly  decrease  the  attainable  null  depths  for 
that  array. 

When  factors  other  than  receiver  frequency  response  do  not  limit  the 
null  depth,  the  general,  more  direct  approach  using  the  FFT  considered  herein 
for  receiver  frequency  response  correction  would  allow  deeper  nulls  to  be 
achieved.  Such  an  approach  is  more  difficult  to  implement,  but  with  state— 
of-the-art  technology  it  can  be  done  for  3  kHz  bandwidths  and  a  four-element 
adaptive  array.  It  should  therefore  be  considered  when  other  limitations  do 
not  already  restrict  the  null  depth,  and  greater  null  depths  are  desired. 
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Writing  the  frequency  response  of  the  receivers  as 
[^(f)  =  rn(f)exp(j4>n(f ))  and  considering  only  random  differential  frequency 
responses  in  the  amplitude  rn(f)  and  phase  4>n(f)  response  of  the  receiver, 
we  get 

/  x  (f)  X 

- « (vf>w„vf>r  v»  (~rm +  3Svf>) 

=  vn(f)  (6pn(f)  +  j64>n(f)) 

where  6pn(f)  is  the  normalized  difference  6r  (f)/rn(f)  in  the  amplitude 
response  from  the  mean  at  that  frequency,  and  <5<j>n(f)  the  difference  in  phase 
response  from  the  mean  at  that  frequency.  The  power  received  in  the  null 
direction  then  becomes 

N 

PN  =  SZ  |Vn(f)|2(62pn(f)  +  62<{.n(f))df. 

n=l 

For  completely  uncorrelated  signals  of  similar  amplitude  in  the  output 
channels,  the  output  power  would  have  been  the  average  array  gain 


P  =  I |v(f) |2df  =  SZ  | V  (f) |2df. 

A  n=l  n 

Using  this  as  a  reference  level,  the  null  depth  is  given  by 


VPA 


n.  |v„(f)|2j«2pn<f)  +  5!*n(f)|df 

n=l _ 

N 

fl  |Vn(f)|2df 
n=l 


°P2  +  °4»2 


where 


fl  |vn(f) |262pn(f)df 

n=l _ 

N 

SZ  |vn(f)|2df 

n=l 


and  o* 


SZ  |  Vn (f )  |  262<J>n(f  )df 

n=l _ 

N 

SZ  |vn(f)|2df 

n=l 


Op  can  be  thought  of  as  the  variance  in  normalized  relative  amplitude 
response  expressed  as  the  square  of  a  dimensionless  ratio  and  0<j,2  as  the 
variance  in  relative  phase  response  <f>»  expressed  in  (radians2),  distributed 
according  to  the  power  lvn(f)|2  in  the  receiver  output  lines. 


